Optimal operating points of oscillators using nonlinear resonators 
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We demonstrate an analytical method for calculating the phase sensitivity of a class of oscillators 
whose phase does not affect the time evolution of the other dynamic variables. We show that 
such oscillators posses the possibility for complete phase noise elimination. We apply the method 
to a feedback oscillator which employs a high Q weakly nonlinear resonator, and provide explicit 
parameter values for which the feedback phase noise is completely eliminated, and others for which 
there is no amplitude-phase noise conversion. We then establish an operational mode of the oscillator 
which optimizes its performance by diminishing the feedback noise in both quadratures, thermal 
noise, and quality factor fluctuations. We also study the spectrum of the oscillator, and provide 
specific results for the case of 1// noise sources. 

PACS numbers: 05.45.-a, 85.85.+J, 62.25.-g, 47.52. +j 
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I. INTRODUCTION 

The importance of oscillators is widely expressed in 
both natural and engineered systems. Such devices, gen- 
erating a periodic signal at an inherent frequency, are 
often developed to serve as highly accurate time or fre- 
quency references [1] . Since the ideal self-sustained oscil- 
lator is mathematically described as a periodic solution 
of a set of autonomous differential equations, its steady 
state dynamics can be expressed in terms of the phase 
variable tracing the motion along a limit cycle. This 
phase is highly sensitive to the stochastic noise present 
in the physical system, as the appearance of periodicity 
without an external time reference implies the freedom to 
drift along the phase direction. The noise induced phase 
drift causes an unwanted broadening of the oscillator's 
spectral peaks, which ideally would have been perfectly 
discrete [2]. Therefore, the phase sensitivity to noise, 
which quantifies the oscillator's performance under the 
influence of an arbitrary noise spectrum, should be as 
small as possible. 

In this paper we study oscillators that are well de- 
scribed by equations of motion that are independent of 
the oscillator phase variable, so that the phase variable 
does not affect the time evolution of the dynamical sys- 
tem. Common examples that fall in this class are oscil- 
lators near the Hopf bifurcation generating the periodic 
motion, and the standard reduced envelope description 
of oscillators based on a high Q linear or weakly nonlin- 
ear resonator, with the oscillations sustained by energy 
injection either through an active feedback loop [2-5], or 
by nondegenerate parametric excitation [6]. An impor- 
tant consequence of this feature is that the limit cycle 
has circular shape, and the complete stochastic problem 
can be resolved analytically as we show in Sec. II. In 
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the absence of this symmetry, the stochastic problem can 
only be treated numerically, using methods such as those 
proposed by Demir [7-10]. Besides this computational 
advantage, a remarkable physical property of this class 
of oscillators is the possibility for complete noise elim- 
ination by tuning a single parameter, which cannot be 
achieved in the general case. 

In Sec. Ill, as an example we calculate the phase sen- 
sitivity to various noise sources of the standard feedback 
oscillator based on a high Q weakly nonlinear resonator. 
It has previously been shown that for a saturated feed- 
back amplitude this oscillator is insensitive to the feed- 
back phase noise at the critical Duffing point [3, 4]. We 
show that this property is a particular case of param- 
eter noise elimination, and provide explicit expressions 
for two feedback loop phase-shift values as a function of 
the amplifier saturation level, for which feedback phase 
noise is completely eliminated, and an additional value 
that eliminates the amplitude-phase noise conversion by 
the nonlinear resonator. We show that in the high satu- 
ration limit, one of the operating points that eliminates 
the feedback phase noise also eliminates noise in the feed- 
back amplitude, additive thermal noise, and quality fac- 
tor fluctuations. Operating at this point at the high sat- 
uration limit is thus an optimal operational state of the 
oscillator, as it was recently demonstrated using an os- 
cillator based on a nanomechanical resonator [11]. 

We also use the simple analysis possible for phase- 
symmetric oscillators to derive complete expressions for 
the spectrum of the oscillator for both white and colored 
noise sources. For the special case of 1// noise sources, 
often encountered in physical implementations, we show 
that the near carrier spectrum is well approximated by a 
Lorentzian for weak noise, and by a Gaussian for stronger 
noise. In both cases the tails of the spectral peaks away 
from the carrier frequency ujq fall off as (uj — ljq)~ 3 , as is 
known from previous results. 



2 



II. NOISE PROJECTION METHOD 

We start by describing an analytical method for calcu- 
lating the stochastic phase sensitivity of a class of oscil- 
lators described by the N equations of motion 



(1) 



where X T = (Xl, Xjy-i, $) is the vector of dynami- 
cal variables, and f T = (/i, /jv)- Note that f is not 
a function of the phase variable $. For the oscillat- 
ing state, the first N — 1 equations have a fixed point 
(Xi.o, -Xjv-i.o), an d the phase continuously rotates at 
the rate wo = /jv(^i,0) •••>-Xjv-i,o)- This is the deter- 
ministic description of a class of oscillators for which the 
phase variable $ does not affect the time evolution of 
the dynamical system. Once we add a particular small 
random noise defined by the stochastic scalar S(i) to this 
deterministic description and linearize the equations near 
the fixed point describing the limit cycle, the dynamics 
of the small perturbations x T = (xi, aijv— i, 4>) are de- 
scribed by the equations 
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witn v noisc(^i,o, — i-Xjv-i,o) = (unoUe,i) "noise, w) de- 
scribing how the noise couples to the dynamical variables, 
and the Jacobian 
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this result is that in the linear approximation in the rotat- 
ing frame, each point on the limit cycle $o is associated 
with a hyperplane H& spanned by the eigenvectors of 
the Jacobian corresponding to negative eigenvalues, and 
all points in this hyperplane asymptote to $o- There- 
fore, any perturbations in this hyperplane have no long 
term effect on the phase. The zero-eigenvalue eigenvec- 
tor of the transposed Jacobian is orthogonal to this hy- 
perplane, and thus projecting the noise vector onto this 
vector accounts for the noise induced phase drift. Under 
the full nonlinear flow, the set 7$ in the stable mani- 
fold of the limit cycle that asymptotes to <£>o is called the 
isochrone [9, 12, 13], and H§ is its linear approximation 
at the vicinity of the limit cycle. 

In Sec. IV we use Eq. (5) and express the oscillator 
spectrum as a function of the noise spectrum. However, 
the phase noise will be proportional to the phase sensitiv- 
ity (v^v no i S0 ) 2 regardless of the particular spectral shape. 
This quantity thus provides an experimentally controlled 
measure of the oscillator performance, under the influ- 
ence of an arbitrary noise spectrum. The only change 
necessary for applying Eq. (5) to a general limit cycle, 
is to replace the two vectors with their time dependent 



(2) counterparts vj_v r 



v ± (w i + 0(t)) J 



(oj t 



-Xl t O,---,-XjV — 1,0 



4>{t)), as derived by many authors [7, 9, 14]. The fact that 
in the present case these vectors are constant leads to a 
particularly simple derivation and analysis, and grants 
the system the possibility of complete noise elimination 
by tuning a single parameter to a point at which the 
scalar product is zero, as we show next. 



Parameter Noise 



Since the oscillator spectral peak is typically much nar- 
rower than the corresponding line width of the driven res- 
onator, which is determined by the decay rates that are 
the negative eigenvalues of the Jacobian, we are mostly 
interested in frequency offsets satisfying \u — u)q\ <C |Aj|. 
This requirement is equivalent to taking the time deriva- 
tives of the first N — 1 variables of Eq. (2) to be zero. If 
we then multiply Eq. (2) from the left by some arbitrary 
vector v T = (i>i, ...Ujy-i, 1), we obtain a scalar equation 
for the time evolution of the phase perturbation 
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Now taking v = vj_, where vj_ is the eigenvector of the 
transposed Jacobian that corresponds to the zero eigen- 
value, for which J T vj_ = 0, eliminates the first term on 
the right hand side of Eq. (4), and we are left with 



,S(t). 



(5) 



Thus we find the result that the noise driven phase evo- 
lution is determined by the scalar product of the noise 
vector and the zero mode of the transposed Jacobian. 
This scalar product quantifies the phase sensitivity to a 
particular noise vector. The geometrical interpretation of 



Let us consider noise which derives from fluctuations 
in some parameter p of the equations so that v^ oise = 
(dfi/dp, dfisr/dp). The fixed points satisfy the fol- 
lowing equations 
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— /jv-i(Xi,0j — , Xn-i,o;p), 
= /jv(-X"i,0i — >-Xjv-i,o;p)- 

First, we find the effect of a change in the parameter p 
on the mean oscillator frequency by differentiating these 
equations and multipling from the left by to find 



d,Lu a = v^JdXo + vXv no i S edp, 



(7) 



with the same vector v no i se appearing as in the stochastic 
problem for noise in the parameter p. Since the first term 
on the right hand side of Eq. (7) is zero we get 



duo 
dp 



' X. noise ■ 



(8) 



Now returning to the stochastic problem, we see that the 
measure which quantifies the phase wandering as a result 
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of parameter noise, is just the derivative of the frequency 
with respect to the parameter containing the noise. This 
implies that for the oscillators under consideration, phase 
fluctuations due to noise in a particular parameter can be 
completely eliminated by tuning the system to the point 
at which dioo/dp = 0. 



Eq. (10) into real and imaginary parts gives the equations 



da 
~dT 

dT 
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sin A = f a (a), 
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(11) 



= U{a). 



B. Amplitude-phase noise conversion 

Other noise sources, such as additive thermal noise, 
will not in general be equivalent to a parameter fluctua- 
tion, and typically it will not be possible to eliminate the 
resulting phase noise completely. To sec how to reduce 
the effects of such a noise source, we write the stochastic 
equation for the phase variable in Eq. (2) explicitly 



JV-l 

E 



(-^1,0, — ) -Xjv-i,o)£. 



(9) 



In correspondence with the general description of the 
previous section, the first equation has the fixed point 
a Q = ssinA, the oscillation frequency is f2 = /*( a o)i 
and the projection vector given by the zero mode of the 
transposed Jacobian is 




cos A 3 
— - 5 _ + - ssi nA,l . (12) 
s sm A 2 



We can now analyze the oscillator phase noise due to 
various anticipated noise sources. 



Feedback noise 



While the last term originates from the noise force acting 
directly on the phase variable, the first N — 1 terms ac- 
count for the phase drift due to noise conversion from 
other fluctuating variables x n to the phase. For os- 
cillators based on a high-Q resonator, amplitude-phase 
noise conversion in a simple two-variable description is 
often depicted as limiting the operating amplitude, since 
this source of phase noise increases when the resonator 
is driven into the nonlinear regime. However, we now 
see that the conversion of noise from the variable x n 
to the phase is completely eliminated by operating at 
duio/dX n = 0. Note that this is a partial derivative, in 
contrast with the full derivative that appears in Eq. (8). 



Noise in the saturated feedback signal is expressed in 
the complex amplitude description by two noise vectors, 
one representing fluctuations in the feedback amplitude, 
and the other representing fluctuations in the feedback 
phase. Since the noise in the feedback phase is more 
dominant for a saturated feedback [3, 4], as a first exam- 
ple of our formalism we calculate the phase sensitivity 
to this noise source, which is expressed by fluctuations 
in the parameter A. The corresponding noise vector is 
= (df a /dA,df$/dA). Using the explicit form of 
va,v_l or the expression derived in §11 A gives 



T 

Vj_V A 



dflo 
~dA 



3s 2 sin A cos A 



1 



2 shr A 



(13) 



III. FEEDBACK OSCILLATOR BASED ON A 
NONLINEAR RESONATOR 



Precise time and frequency references are typically 
built around high Q resonators, since high-Q reduces 
phase noise in the linear regime [15]. In the high Q limit, 
the standard feedback oscillator is effectively described 
by a complex amplitude equation that captures the slow 
dynamics near the high natural frequency u)q of its weakly 
nonlinear resonating element [2, 3, 5, 16] 



dA ( 1 ■ 3 ,,,2 

dT= -2 + *8 |A| 
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where T = tuio/Q is a slow time that scale, A = 
the complex amplitude, A is the phase shift of the feed- 
back, and s is the feedback level, which is assumed to 
be the output of a saturated amplifier and so is indepen- 
dent of the magnitude |A| of oscillation [25]. Separating 



At two phase shift values A 1,2, this expression is zero 
and the oscillator is insensitive to the phase noise of the 
feedback signal. The explicit values of Ai^ are given in 
Eq. (A7) in the Appendix. These points bifurcate at the 
critical Duffing point for which Ai = A2 = 27r/3. As we 
already noted, the ability to eliminate the feedback phase 
noise at the critical Duffing point has been previously 
established [3, 4]. 

The second contribution of the feedback noise to the 
phase drift comes from noise in the feedback level. This 
noise is expressed by fluctuations in the parameter s, for 
which the noise sensitivity is 



ds 



3s sin A 



(14) 



This noise cannot be eliminated in the same manner as 
feedback phase noise, since the oscillator amplitude is 
zero at A = rm, however we show in §111 C how it can be 
significantly reduced. 
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two quadratures. In the complex amplitude representa- 
tion, this is described by noise vectors for the two sources, 
= (1,0), and = (0,1/a). The total phase sensi- 
tivity to thermal noise is given by adding the two noise 
terms in quadrature 



TPS(A) = (v^vcj) 2 + (vlv K ) 5 



(18) 



1 



s 2 sin A 



6s . 
— sin A 

4 



cos A 
s sin 2 A 



Since the hrst term is the direct phase noise, and the sec- 
ond results from amplitude-phase conversion, the second 
term satisfies v^Vjf(A a |0) = 0. However, the total ef- 
fect of thermal noise is minimal at the nearby phase shift 
value Atps, which satisfies rfTPS(ApsT)/^A = 0. 



FIG. 1: (Color online) The amplitude vs. frequency curve of 
the oscillator described by Eq. (10). The low and high am- 
plitude red dots are Ai,2, respectively, for which dQo/dA — 
dflo/da — 0. The star is A^ for which Otto /da — 0. The 
bottom curve is for s = 3, and the top one is for s = 10, in- 
dicating that Ai and A a |^ approach each other as saturation 
level increases. 



B. Amplitude-phase noise conversion 

Based on the analysis of §11 B, complete elimination 
of amplitude-phase noise conversion is achieved in the 
present oscillator for 



da 



scosA 3s 1 cos A n ,„_. 

- r = -BinA+- ,. =0, 15 
ai 4 2 S sin 2 A 



whose real solution A a |0 is given in Eq. (All) in the 
Appendix. As the saturation level grows AgU approaches 
Ai as illustrated in Fig. (1), which shows these points on 
the amplitude-frequency curve for the oscillator, which 
follows the well known Duffing curve. Note also that 
from Eq. (13), we have 



gOo(Ai) 
da 



1 



2s cos(Ai) ' 



(16) 



so that as the saturation grows, operating at Ai elimi- 
nates the amplitude-phase noise conversion as well. On 
the other hand, at the operating point A^a we have 



drj (A £ 



dA 



1 

2' 



(17) 



and the feedback phase noise is not reduced here. Thus, 
although the points A a |0, Ai approach each other as sat- 
uration grows, it is better to operate at Ai than at A a |^. 



C. Optimal operating point for combined noise 



sources 



The four special phase shift values minimizing the ef- 
fects of various noise sources on the feedback oscillator 
are plotted in Fig. 2, as a function of the other experimen- 
tally controlled parameter, the saturation level. Since the 
real physical oscillator is subjected to multiple, uncorre- 
cted noise sources, we would like to use these results to 
find its optimal operational state. We know from §111 A 
that the phase shift values Ai^ eliminate the feedback 
phase noise. Since this is a substantial noise source we 
would like to eliminate, lets examine the sensitivity to 
feedback amplitude noise at these points. Using Eq. (14) 
and the high saturation limit of Eq. (A6) allows us to 
deduce the limits 



lim v^v s (A 2 ) ~ s ->• oo, 
lim v^v s (Ai) - s- 1/3 -> 0. 



(19) 



Therefore, Ai is a preferable operating point from the 
feedback noise perspective because for high saturation 
it eliminates noise in both quadratures of the feedback, 
while operating at A 2 only eliminates the feedback phase 
noise, with the feedback amplitude noise becoming more 
important with increasing saturation level. 

To consider the effects of thermal noise at Ai, we 
use Eq. (16), which implies that for large saturation 
amplitude-phase conversion is diminished at this point, 
and together with the limit lim^oo ao(Ai) ~ s 1 / 3 , we 
derive the results 



lim v[vr(Ai) ~ s- 1 -> 0, 



(20) 



lim v^vq(Ai) ~ s 



'1/3 



0. 



1. Thermal noise 

Thermal noise acting on the resonator degrees of free- 
dom is manifested as equal independent sources in the 



Thus, in the large saturation limit, operating at Ai also 
eliminates both components of thermal noise. 

Noise in the quality factor is expressed by fluctuations 
in the linear damping coefficient in Eq. (10). Since we 
consider the damping to be small, fluctuations with the 
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FIG. 2: Special phase shift values as a function of the 
drive amplitude. Ai,2, which eliminate the feedback phase 
noise, bifurcate at the critical Duffing drive amplitude s = 
(4/3) (5/4) . On increasing saturation level, A2 approaches tt/2, 
and the other three special phase shift values approach n. 



same relative intensity as other parameters will typically 
yield a considerably smaller effect on the phase drift. 
However, it is instructive to include these fluctuations 
in the analysis as well. The noise vector representing 
them is Vq = (—a/2,0); it is transferred to phase noise 
through amplitude-phase conversion and satisfies 

vlv Q (A a | ) = 0, (21) 
lim v;[v Q (Ai) - s~ 2/3 -> 0. 

Furthermore, since generally damping terms do not ap- 
pear in the dynamic equation for the phase variable, fluc- 
tuations in nonlinear damping coefficients arc also elimi- 
nated at A a |0. For the often encountered case of nonlin- 
ear damping proportional to the amplitude cubed [16- 
18], the phase sensitivity at the high saturation limit 
scales as a constant (~ s°) at Ai. 

Therefore, we have shown that increasing the satura- 
tion level of the amplifier and setting the phase-shift to 
Ai, which in this case approaches 7r, optimizes the os- 
cillator's performance. This operational mode eliminates 
the contributions of the feedback noise, thermal noise, 
and quality factor fluctuations on the oscillator's phase 
drift. The phase sensitivity to combined noise sources, 
whose minimum is obtained for a phase shift value that 
approaches Ai as saturation grows, and the breakdown 
to the contributions of the different noise sources, are 
shown in Fig. 3. These results were recently verified by 
experimental phase noise measurements of an oscillator 
based on a nanomcchanical resonator [11]. 

Additional sources of phase noise might originate from 
the elasticity mechanism of the resonator [19]. These are 
expressed by fluctuations in the resonance frequency and 
nonlincarity coefficient. Fluctuations in the resonance 
frequency directly translate into phase noise, whose level 




FIG. 3: (Color online) Phase sensitivity as a function of the 
phase shift, (a) Different contributions (v^v no i sc ) 2 to the 
phase noise, for s = 3. We see that the feedback phase noise 
curve is zero at Ai,2, and that those of the real component of 
thermal noise and quality factor noise are zero at A a |^, which 
is slightly to the right of Ai . (b) The total phase sensitivity, 
which is a sum of all the curves in (a), for different values 
of the drive level. We see that as saturation level grows, the 
optimal operating point approaches Ai, noted as red asterisk, 
and the noise level at this point decreases. 

is approximately constant in the experimental parame- 
ters. Fluctuations in the nonlinearity coefficient grow as 
a function of the saturation level; however, if the fluctu- 
ations in both these parameters are correlated, we can 
potentially design the resonator nonlincarity in such a 
way so that these contributions cancel each other out. 



IV. THE OSCILLATOR SPECTRUM 

The simplicity of Eq. (5) allows us to derive the full 
oscillator noise spectrum for a variety of noise sources 
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(white, colored, 1/f. . . ). Many of these results may be 
found in the previous literature [1, 2, 7, 8, 20, 21], but the 
present approach gives a particularly compact and gen- 
eral formulation for the class of oscillators we consider. 
We assume a stationary Gaussian noise source with the 
autocorrelation (2(£ 2 )S(fi)) = Rjyfo—ti), and the spec- 
trum given by a fourier transform Sn(w) = •7 r [-Riv(i)]- 
Then <f){t + t) - 4>{t) = P J* +T E{t')dt' (r > 0) is also 
Gaussian with P = VjV no i se the noise projection con- 
stant, and variance V(t) = ((<f>(t + t) — <f>(t)) 2 ) given by 



V(r) = 2P 2 / (r-v)R N {v)dv 
Jo 



4P 2 

— / 



sin(a;r/2) 



du>. 



(22) 



The variance V(r) typically grows with r: for a noise 
source with correlation time r c [26], for large enough 
r 3> t c the variance increases linearly with time V(r) ~ 
P 2 Sn(0)t corresponding to phase diffusion. 

The spectrum of the oscillator is the Fourier transform 
of the autocorrelation function (D(t + r)D(t)), where 
D(t) is the displacement of the resonator at time t. We 
take this to be given by D{t) = acos(u)ot + 0). Neglect- 
ing amplitude fluctuations and after transients have died 
out, the oscillator correlation function is [2, 7, 22] 

(D(t + r)D{t)) ~ ^! e -^)/ 2 cos(w r) = R{t) (23) 

where the Gaussian properties of <fi and the identity 
[e lx ) — e~' x '' 2 for a zero mean Gaussian stochastic 
variable x have been used. We see that this autocor- 
relation function of the oscillator is indeed a station- 
ary stochastic process, as expected for a system with- 
out an external time reference. The oscillator noise 
spectrum S(uj) = J-[R(t)] can be written as S(u>) = 
a 2 [S(uj + uj ) + S{uj - w )]/4 with 



S(u) = T\e 



-V(t)/2 ] 



(24) 



The oscillator spectrum about the carrier frequency is 
thus given by evaluating this Fourier transform. 

For large frequency offsets, the Fourier transform in 
Eq. (24) is dominated by small times where the variance 
V(t) is small. In this case the exponential can be ex- 
panded to first order, so that for w ^0 



\F[V(r)]=R 



S N {u>) 



(25) 



This corresponds to the well known Leeson results [15] 
for the oscillator noise spectrum, namely an lo~ 2 depen- 
dence for a white noise source, w~ 3 for 1/f noise etc. On 
the other hand, for u — > 0, the Fourier transform is de- 
termined by the long time behavior when the variance 
V(t) gets large, and so the full exponential expression 
must be used. For small enough noise, V(t) ~ 1 only 



at long enough times such that the diffusion behavior 
V(r) oc r applies. In this case the Fourier transform 
gives a Lorentzian, so that for these frequencies 



4P 2 S N (0) 



4u 2 + [P 2 S N (0)Y 



(26) 



The frequency where the spectrum crosses over from the 
Lorentzian Eq. (26) to the power law tail Eq. (25) will be 
of order f _1 where f is the time at which the variance 
V(t) grows to be 0(1). For a white noise source the 
Lorentzian leads directly to Eq. (25) at large frequencies. 

Oscillator phase noise for a general limit cycle and 
white or colored noise has been looked at previously. For 
a general limit cycle both the noise vector and projection 
vector will be time dependent around the limit cycle. Al- 
though the oscillator spectrum in the general case is still 
given by (24), the corresponding phase variance changes 
to 



(27) 



Vgcnoral(T) = / dll dvR N (v) 



x v^(w)v noisc (u)v^(w - v)v noisc (u - v). 

This much more complicated expression (cf. Eq. (22)) 
has only been evaluated in certain limits. Demir et 
al. [7] looked at the case of white noise (r c — > 0) for 
which expressions Eqs. (25,26) apply with the replace- 
ment P 2 -> t~ x / ' p (v^(i)v noise (t)) 2 G?t (ie. the mean 
square value of the noise projection). On the other 
hand Dcmir [8] looked at the case of colored noise 
with correlation time much longer than the oscillator 
period r c 3> t p , deriving the results Eqs. (25,26) but 
now with the replacement P — s- f* p v^(i)v no i se (£)(i£ 
(ie. P 2 is replace by the square of the mean of the 
noise projection around the limit cycle). Nakao et 
al. [14] look at general colored noise, but only derive 
results for the Lorentzian component of the spectrum 
(not the tails further away from the carrier) given by 
the long time, diffusive behavior of V(t). They find 
a more general expression for the diffusion constant, 
leading to the replacement in Eq. (26) P 2 Sn(0) — > 

dvR N (s) t' 1 J Q P duv^(w)v noiso (u)v[[(« 
v)v no j se (w — v). For the special class of oscillators 
considered in the present work, both v no ; sc ,v^ arc 
constant and our expression (22) encompasses these 
limiting cases. 



A. Oscillator spectrum for 1/f noise 

An interesting special case that arises in may physical 
implementations is a 1/f noise source (or more generally, 
l/f(a-i) w i t h 1 < a < 2). Using the integral 



o (2tt/)« + x a 



dx = 



a sin i 



(/" 



(28) 
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a convenient representation of a l//^" 1 ' spectrum, with 
a low frequency cutoff to eliminate the divergence as / — > 
and give a finite correlation time is 



Sn(I) oc 



1 



(2tt/) q + x a 



dx. 



(29) 



We focus on the pure 1/f case given by a = 2, so that 
the noise spectrum is 



S N (f) = 4/ 



dx 



x 2 + {2nf) 2 
I 2 4/ 2 arctan(w c /(27r/)) 

~ m 2^7 ■ 

The corresponding phase variance is [8] 
2P 2 I 2 - 



(30) 



VM 



+ (tuj c ) 2 E(tlu c ) 



(31) 



with E(x) = J^ 00 dye xy /y the exponential integral func- 
tion. Note that for r ^> a;" 1 , V(t) oc t giving phase 
diffusion, but for shorter times r < lu// 1 



V(t) ckt 2 (I - 7 - Iog(w c r)) , 



(32) 



where 7 ~ 0.58 is the Euler-Maschcroni constant. 

For very weak noise, the arguments of the preceding 
section apply, leading to the results Eqs. (26,25), ie. a 
Lorentzian with u~ 3 tails. However since we might ex- 
pect the cutoff oj c to be small, even for not too large noise 
strengths it is possible for the phase variance V(r) to be- 
come comparable to unity for times still small enough 
for Eq. (32) to apply rather than the diffusive behavior. 
The resulting Fourier transform Eq. (24) can be approx- 
imately evaluated in this case by replacing the log(w c r) 
term in Eq. (32) by the constant log(w c /oj m ) where u> m 
represents the frequency range over which we want to ap- 
proximate. The corresponding spectrum is a Gaussian 



(33) 



with a 2 = P 2 I 2 (3/2 - 7 - \og{w c /w m )). The resulting 
oscillator spectra for weak and strong noise intensities are 
shown in Fig. 4. A Gaussian spectrum near the carrier 
frequency also results from the noise spectrum 1//' 1_£ ^ 
with < e -C 1 [23], and for a rectangle noise spectrum, 
a Lorentzian and Gaussian are accepted in the weak and 
strong noise limits, respectively [24]. 



V. CONCLUSIONS 

We have demonstrated an analytical method for cal- 
culating the phase sensitivity of a class of oscillators, 
which allows for complete elimination of certain phase 
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(a) 



(b) 



FIG. 4: The spectrum of an oscillator subjected to 1/f noise. 
The exact solution is T[e~ v{T)/2 ] with V(r) given by Eq. (31). 
For very weak noise, small offset frequencies are described by 
the Lorentzian (26), however for stronger noise, this frequency 
regime is approximated by the Gaussian (33) as shown in (a) 
and (b), respectively. Both curves approach the expression 
(25) at large frequencies. Other parameters are uj c = 0.1, 
P 2 = 1, and uj m = 10. 



noise attributes. We have applied this method to the 
standard feedback Duffing oscillator, and exposed an op- 
erational mode of this oscillator which diminishes both 
components of the feedback noise and thermal noise, and 
noise in the damping mechanisms. This therefore offers 
an optimal operational state of the oscillator which uses 
the resonator nonlinearity to reduce its phase noise, and 
disproves the common perception that amplitude-phase 
noise conversion necessarily degrades the oscillators per- 
formance at the nonlinear regime. We have also studied 
the full oscillator noise spectrum for white and colored 
noise sources. For a white noise source, or noise with 
a short correlation time, the near carrier spectrum is 
Lorentzian, and so varies as (w — loq)~' 2 towards larger 
frequencies. For the often encountered case oil/ f noise 
the near-carrier spectrum is a Lorentzian for weak noise, 
but as the noise level grows it is well approximated by a 
Gaussian; in both cases the phase noise spectrum crosses 
over to (to — luq)~ 3 away from the carrier frequency. 

This research was supported by DARPA through the 
DEFYS program. 



Appendix A: Explicit expressions for special phase 
shift values 



1. Feedback phase noise 

The phase shift values for which feedback phase noise 
is eliminated are given by equating Eq. (13) to zero, 



3s 2 sin A cos A 1 



2 sin 2 A 



0, 



which is equivalent to 



x 4 - 2x 3 + 2x + — - 1 = 0, 



(Al) 



(A2) 



8 



with x = cos 2A, and C = 9s 4 /16. The two real solutions 
of this equation are 



X2 



(l-r + «), 
(l-r-q), 



(A3) 



with r = y'l + y R , q = 
real root of the equation 



V 



whose discriminant is 



y/Z — r 2 + 2/r, and yn being a 

(A4) 



(l + t/)^ = 0, 



D = - 



256(27C- 64) 
C3 ' 



(A5) 



Since the coefficients of Eq. (A4) are real, for < C < 
64/27 it has three real roots; however two of them are 
negative and smaller than -1, so r is imaginary. For 
the third, positive root, on increasing C the value of q 
changes from imaginary to real at C = 64/27, which is 
the critical Duffing point, corresponding to the drive am- 
plitude s = (4/3) 5 / 4 [16]. The corresponding root is 



VR 



2 3 1 / 3 ^3C 3 (27C - 64) + 9C 2 + 3 2 / 3 4C 



2/3 



'3C 



3C ( v / 3C 3 (27C-64) + 9C 2 
cosh 



1/3 



^rctanhf ^ 3C3(2?C ' 64) 
3 arctanh ^ 



(A6) 



and the phase shift values which solve Eq. (Al) are 

arccos(xi 2) 

A li2 =7T 

2. Amplitude-phase noise conversion 



(A7) 



Equation (15) is equivalent to 

s o 1 + x 

(1 - xf - — = 0, 

whose negative discriminant is 

4(27C + 1) 



(A8) 



V = -- 



c 3 



(A9) 



Equation (A8) thus has one real solution, which is 



( v /3C 3 (27C+l)-9C 2 ) 1 / 3 

1 + 3 2 / 3 C 
1 



V / 3(v / 3C 3 (27C+l) - 9C 2 ) 1 / 
and Eq. (15) is thereby solved by 

arccos(a; a |0) 



T73 ' ( A1 °) 



A4> 



(All) 
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